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Abstract 

We consider a system of N nonrelativistic particles of spin 1/2 
interacting with the quantized Maxwell field (mass zero and spin one) 
in the limit when the particles have a small velocity. 

Two ways to implement the limit are considered: c — > cxd with the 
velocity v of the particles fixed, the case for which rigorous results 
have already been discussed in the literature, and u — > with c fixed. 
The second case can be rephrased as the limit of heavy particles, 
mj s~'^mj, observed over a long time, t £~^t, e O"*", with 
kinetic energy E'kin = 0(1). 

Focusing on the second approach we construct subspaces which 
are invariant for the dynamics up to terms of order e-y/log(e~^) and 
describe effective dynamics, for the particles only, inside them. At the 
lowest order the particles interact through Coulomb potentials. At the 
second one, e^, the mass gets a correction of electromagnetic origin 
and a velocity dependent interaction, the Darwin term, appears. 

Moreover, we calculate the radiated piece of the wave function, i. 
e., the piece which leaks out of the almost invariant subspaces and 
calculate the corresponding radiated energy. 

1 Introduction 

A system of nonrelativistic particles of spin 1/2 interacting with the quan- 
tized radiation field is described by the so-called Pauli-Fierz Hamiltonian, or 
"nonrelativistic quantum electrodynamics". The model is thought to have 
an extremely wide range of validity, apart from phenomena connected to 
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gravitational forces and from other ones typical of high-energy physics like 
pair creation, whose description requires the use of full relativistic QED. 

This belief is mainly based on the analysis of some formal limit cases, 
which can be accurately studied both from a theoretical and an experimen- 
tal point of view. Indeed, the interaction between charged particles is usually 
described by instantaneous pair potentials of Coulomb-type, without intro- 
ducing the field as dynamical variable. This is known to be a good approx- 
imation if the particles move sufficiently slowly. One aim of this paper is 
a mathematically rigorous justification of this fact, i. e., the derivation of 
the Schrodingcr equation with Coulomb potentials, and the second order ve- 
locity dependent corrections to them, starting from nonrelativistic quantum 
electrodynamics. In addition, a formula is provided for the wave function 
of the radiated photons and the corresponding radiated energy, which is the 
quantum equivalent of the Larmor formula of classical electrodynamics. 

In more detail the model considered is given, excluding the addition of 
the electronic spin, by the canonical quantization of a system of N classical 
charges interacting through the Maxwell field. 

A sharp ultraviolet cutoff is introduced assuming that each charge has a 
charge distribution given by 

Qj{x) = ejip{x), X E W^, (1) 

where the form factor satisfies >^{k) — (27r)~^/^ for < A, otherwise (note 
that there is no infrared cutoff'). 

The classical equations of motion are given by 

-d^B{x,t) = -V X E{x,t), 



c 



1 ' (t) 

-d^Eix, = V X B{x, t)-J2 ^JV{^ - ^ 



ejcp{x-qj{t))'^, 



with the constraints 



N 

V • E{x, t)^Yl ^^^(^ - W) ' ^ ■ ^(^' = 0' (3) 
and the Newton equations for the particles, 

miqiit) = ei[E^{qi{t),t)) + ^ X B^{qi{t),t)], l^l,...N, (4) 



where E^{x,t) := {E *x and analogously for 5, 
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The canonical quantization of this system in the Coulomb gauge is de- 
scribed, e. g., in ( |Sp| , chapter 13). The Hilbert space of the pure states is 
given by 

,^:=^^®^. (5) 
The space for the particles, is defined bjfl 

:= L2(M3 X Zs)®"^, (6) 

where is the configuration space of a single particle and Z2 represents its 
spin. 

The state space for a single photon is //^(R'^ x ^2), where M'^ is the mo- 
mentum space of the photon and Z2 represents its two independent physical 
helicities. The photon Fock space is therefore 

=^:=©M=o<)^'(K'x^2), (7) 

where denotes the M-symmetric tensor product and x Z2) := 

C We denote by i^p the vector (1, 0, . . .), called the Fock vacuum. 
The dynamics of the system are generated by the Hamiltonian 

2 

+ V"^coui(a;) + ciyf, (8) 

where all the operators appearing are independent of c and we use units in 
which ^ = 1. (Tj is a vector whose components are the Pauli matrices of 
the jth particle, A^p{xj) denotes the quantized transverse vector potential 
in the Coulomb gauge, V^coui is the smeared Coulomb potential and iff the 
free field Hamiltonian. The reader who is not familiar with the notation is 
advised to look at section [2], where the model is described in more detail. 

To implement practically the idea that the particles move "slowly", a 
standard procedure, applied also in classical electrodynamics (see, e.g., pa] . 
|LaLij ). is to take the limit c ocEl. Since c is a quantity with a dimension, 
one should actually say that |f |/c — > 0, where f is a typical velocity of the 

^The formalism presented holds also in the case when all the particles are equal and 
their Hilbert space is given by the subspace of totally antisymmetric wave functions. In 
this case, the dipole radiation given in (j29p is zero. We consider therefore the general case 
of different particles. 

^In the classical case, a more refined and precise analysis is carried through in |KuSpi| , 
|KuSp2| . The authors consider, loosely speaking, initial conditions which represent free 
particles moving together with the field they generate ( "dressed" particles or charge soli- 
tons), with a velocity of order 0[e^I'^') with respect to the speed of light. Assuming that 
the particles are at time t = far apart (relative distance of order ©(e"^)) and rescaling 



^ Orr 



1 2mj 



iV - -j.ejA^{xj] 
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particles. This can be achieved in two ways, fixing v and letting c ^ oo or 
fixing c and letting v —>■ 0. 

In the classical case this is refiected in the fact that the limit c ^ oo is 
equivalent, up to a rescaling of time, to the limit of heavy particles, as one 
can easily verify replacing in equations ©-([I]) rrii with e~'^mi, t with e~^t, 
and looking at the limit e ^ 0. 

We will show that in the quantum case the two procedures are non equiv- 
alent anymore, a fact that can be intuitively explained by the presence of an 
additional scale given by h. 

In this paper we concentrate on the limit of heavy particles observed over 
a long time. An additional aim is to point out similarities and differences 
between the two limits in the quantum context and to compare the results 
we get for the Pauli-Fierz model with the ones valid for the Nelson model, 
where the particles and the photons are spinless, [TeTe] . 

We recall briefiy in the next subsection some known results about the 
case c — cx) and then illustrate in more detail the limit e ^ 0. 



1.1 The limit c ^ oo 

This case, as observed by Spohn |Spj , has the form of a weak coupling limit, 
fact which was already noted for the Nelson model by Davies [Dag], who 



formulated also a general scheme to analyze the limit dynamics in the weak 



coupling case Dai| (an extended notion of weak coupling limit for Pauli-Fierz 



systems in which the Hilbert space of the particles is finite dimensional has 
been examined in [DePe] ) . 

Davies looks loosely speaking at the limit A — > for the time evolution 
generated by an Hamiltonian of the form 

which corresponds physically to a weak interaction, whose effect is however 
observed over the long time scale A~^. 



suitably the dynamical variables, they show that the particles remain at a relative distance 
of order 0{£~^) for long times (of order Oie'^/"^)) and on this time scale their motion is 
governed by effective dynamics. The possibility to implement an analogous limit in the 
quantum case is unclear, because there is no obvious quantum counterpart to the classical 
charge solitons. The Pauli-Fierz Hamiltonian without infrared cutoff has indeed no ground 
state in the Fock space for fixed total momentum different from zero [E] [Chj . We stick 
therefore to the more pragmatic choice c ~* oo. 

The states which we define through the dressing transformation in equation (P7)) 
should be considered approximate dressed states valid for small velocities of the particles. 
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The Hamiltonian assumes a similar form if we consider it on the long 
time scale defined by c^. Putting 

A := 

we get indeed 

c'H' = X-\Hi + X^/^H^ + A/ii + A^/=^/i4/3) := A-^i/;,, (9) 



where 



^ 1 



P ^ 2m, 

j=i 1 

N 

N g2 

V3:=E^^^^(^^)'^ ' (12) 

.7 = 1 

where := V x and we normal order the quadratic term. 

Applying Davies scheme one gets in the end ([Spj, theorem 20.5) 

Theorem 1. Let if) e H^R^^X^""), then 

lim II (e-'^''^'* - e-'^<^— "'*)^/' ® fipllif = 0, (13) 



where 



Hfisrw ■ Hp -\- C V^arw ~l~ C V^pin; (1^) 



darw 



:= - E ^ / ^T^^^'^'^^i^^. ■ (1 - ® K)iV..e-'=--,(15) 
f-^ mjmi 2|/c|2 



N 



where {k ® K)ij := KjKj and k := A;/|A;|. 

Vdarw gives rise to a correction of electromagnetic origin to the mass of 
the particles and to a velocity dependent potential, called the Darwin term. 
It already appears in classical electrodynamics when the dynamics of the 
particles are expanded up to terms of order (f/c)^ (see, e. g., [Jaj or |LaLij ). 
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For the convenience of the reader and to ease the comparison with the 
results for the hmit of heavy masses we give a formal derivation of this 
theorem in appendix |Al 

We note here that the method employed in the weak coupling case forces 
one to consider as initial condition for the field just the Fock vacuum, which 
contains no photons at all. There is therefore no analogy with the physical 
picture that every particle should be described by a "dressed state" , loosely 
speaking the particle itself dragging with it a cloud of "virtual" photons. 



1.2 The limit m ^ oo 

The situation is different in the case rrij oo, which is more conveniently 
studied adopting units where c = 1. Replacing nij by e~'^mj we get then the 
Hamiltonian 

^1 

:= V - — + y^coui + Hi- e-^pj ■ A^{xj) - e^^aj ■ B^{xj) + 
2mi ■> ^ rrii 2m,- 

j — -^ J J J 

(17) 

where we have indicated with pj the e-momentum of the jth particle 

pj := -ieV,^.. (18) 

As already pointed out talking about the classical case, the dynamics 
have to be observed over times of order 0{s^^). This is necessary in order to 
see non trivial effects, because we consider initial states with bounded kinetic 
energy. Since the particles have a mass of order (^(e"^) this means that their 
velocity is in the original time scale of order 0{e). 

To analyze the limit e — > we construct a unitary dressing transformation 
%'s : —>■ J^, which allows us to define dressed states for small velocities 
of the particles and to introduce a clear notion of real and virtual photons. 
More precisely, in the new representation defined by '^e, the vacuum sector 
® flp corresponds to states of dressed particles without real photons, 
while in the original Hilbert space a state with M real photons is a linear 
combination of states of the form 

® a(/i)* ■ ■ ■ a(/M)*fiF), where 7/> G JTp, /i, . . . , /m G L\R' x Z^). 

The projector on the subspace corresponding to dressed particles with M 
real photons is therefore 

Plj:=%*ilp0QM)%, 
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where Qm denotes the projector on the M-particles subspace of the Fock 
space. 

In short, we will show that the subspaces Pm-^ are approximately in- 
variant for the dynamics defined by on times of order 0{e~^). Moreover, 
on this time scale we will give effective dynamics for states inside such a sub- 
space, with an error of order 0{e^\og{e~^)). The effective dynamics contain 
the Darwin correction described in (fTHI) . but no spin dependent term. One 
can get an idea of why this happens comparing the expression of with 
that of Hx, equation ( JTTl) and ([9]). In the spin dependent term is of second 
order, while in H\ is of the first one. In the limit £ — *• the analogue of V^pin 
would be of order therefore it does not appear in an expansion of the 

time evolution till second order. Finally we compute the leading order part 
of the state which makes a transition between and , which corresponds 
to the emission of one real photon. The corresponding radiated energy is 
given by a quantum analogue of the Larmor formula. 

The procedure to construct the unitary is explained in detail in [TeTej 
for the Nelson model. The technique used is based on space-adiabatic pertur- 
bation theory jTej, a method which allows to expand the dynamics generated 
by a pseudo differential operator with an e-dependent semiclassical symbol. 

The main difficulty in all models concerning the interaction of particles 
with a quantized field of zero mass is that, because of soft photons, the 
principal symbol of the Hamiltonian has no spectral gap, which is a condition 
required to apply the methods of |Tej . In the case of we have indeed 



For every fixed (p, q) this is an operator on and has a ground state given 
by fip, at the threshold of the continuous spectrum. The corresponding 
eigenvalue 



is the symbol of an Hamiltonian acting just on and describing the particles 
interacting through the smeared Coulomb potential. 

The trouble connected to the absence of the spectral gap is solved by 
introducing an effective gap, considering the Hamiltonian if^'*^ where the 
form factor (p (see equation ([T])) is replaced by (^cr(^) := (27r)~^/^ for a < 
|/c| < A, otherwise. 

Proposition, (see proposition\3\) 




■p] + K'coui(g) + Hi, (p, q) e M=^^ X M' 



(19) 




(20) 
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Suppose that the cutoff a is a function of e, a = o-{e), such that a{e) < e^, 
then 

ll^-itHVe _ e-^'^-'^/^cm,^.^ < Cmey/' (21) 

where 

J^o-=D{H',) (22) 
is the domain of the free Hamiltonian 

^ 1 

Ho--J:^P^ + H, (23) 

j=l 

with the corresponding graph norm. 

Fixing 0", e. g., as a sufficiently high power of e we can then replace the 
original dynamics with infrared cutoff ones. 

For if^'*^ it is possible to build a dressing operator ^£,a which can be 
expanded in a series of powers of e with cx-dependent coefficients which are 
at most logarithmically divergent. Using it we define the dressed Hamiltonian 

^d;:s ■■= %,^H^'''%:. (24) 

which can be expanded in a series of powers of e in £(t^, ^), with coeffi- 
cients which are also at most logarithmically divergent in a. The different 
coefficients in the expansion correspond to different physical effects which 
can be now clearly separated according to their order of magnitude in e. 

The first result we find, as we already mentioned above, is that the dressed 
M-photons subspaces are approximately invariant for the dynamics: 

Theorem, (see corollary IE)- 

Given a x ^ C'^(]R) and a function a{e) such that 

e-V(e)i/2^0, eVlog(a(e)-i) ^ 0, £ ^ 0+, (25) 

then 

\\[e-''''^,PlMHn\\ci.^) = 0{^^MTT\t\e^/\og{a{ey^)) , (26) 

where 

P!,:=%:,(4lp^QM)%Me) ■ (27) 

The adiabatic decoupling which guarantees the invariance of the sub- 
spaces holds uniformly only on states in which the particles have a uniformly 
bounded kinetic energy. For this reason we introduce a cutoff function on the 
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total energy x, which gives rise automatically to a bounded kinetic energy 
for the slow particles. 

In the following we assume that the function a{e) has been fixed so that 
( l25l) is satisfied. One can then approximate the dynamics of the particles 
inside each almost invariant subspace. 

Theorem, (see theoremlE)- 

Let S be a bounded observable for the particles, S G C{Jifp), and u G 
^i{Pl^x{H^)'^) a density matrix for a mixed dressed state with M free 
photons, whose time evolution is defined by 

u{t) := e-'*^^/^a;e'*^^/^ . 

We have then 

Tr^ [{S ® l^)cj(t)^ = Tr^p (^^e-'*^-ffTr^^(cu)e'*^-ff^ + 

+ 0{e"^\t\){l - 6mo) + 0{6'\og{aie)-'m + \t\')), 
where Smo = 1; when M = 0, otherwise, and 
^ 1 

[ dk ^-^e^'-^p, ■ (1 - ® «:)p.e-^- = (28) 
f-^ rrijmi J ^3 2\k\^ 



i,j=i 1 ' 
^ 1 



2m, 

Remark 1. Even though the subspaces depend on the choice of the in- 
frared cutoff, the effective Hamiltonian is infrared regular and therefore in- 
dependent of a. Moreover, as we briefly mentioned above, it contains the 
corrections to the mass of the particles and the Darwin term, but no spin 
dependent term (compare with theorem [T]). This topic is further discussed 
in the proof of theorem [5] and in remark HI 

Since the subspaces are only approximately invariant, there is a piece 
of the wave function which "leaks out" in the orthogonal complement. This 
correspond physically to the emission or absorption of free photons. For a 
system starting in the dressed vacuum the leading order of the wave function 
of the emitted photon is given in the next theorem. 
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Theorem, (see corollaryl^. 

Up to terms of order 0[e^ log(cr(£:)^"'^)(|t| + the radiated piece for a 

system starting in the dressed vacuum (M = 0) is given by 



where e\{k) is the polarization vector of a photon with helicity X, 

^ 1 

ijix) :=< fip, x(/^o)^ >^ e ^p, (31) 

N 

D{s;x,p) ■.= y^^xf{s;x,p), (32) 
i=i ^ 

Op^ denotes the Weyl quantization acting on a suitable symbol space on 
X and x^' zs i/ie solution to the classical equations of motion 

mjxf{s;x,p) = -Vx,V^coni{x''\s;x,p)), 

xf{0;x,p) = Xj, xf{0;x,p) =pjmj\ j = l,...,iV. 

Remark 2. As explained in detail in remark 0, generically the norm of the 
radiated piece is bounded below by 0[elog{€a{e)~^)), which means that 
the subspace Pq near optimal, i. e. the transitions are at least of order 
0{e\ogiEaiE)-')). 

Note that, like in classical electrodynamics, when all the particles are 
equal, the leading order of the radiated piece vanishes, because D is then 
proportional to the position of the center of mass, whose acceleration is zero. 

Remark 3. Even though the radiated wave function has no limit when e —>■ 0, 
because ip{k)\k\~'^^'^ ^ L^(M'^), the corresponding radiated energy has a limit. 
Defining 

^rad(t) := (^rad(t),^f^rad(t)), (33) 

we get to the leading order (see remark E]) 

^rad(t):=|i?.ad(t) = |^(v,opr(iz)(t)r)^)^, • (2^) 
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In the case of the Nelson model analogous results are proved in |TeTej . 
which contains also a detailed discussion of the adiabatic framework. The 
form of the effective dynamics is equal, the only difference, as one can expect, 
is in the radiated piece, which contains here explicitly the helicity of the 
photon. Another difference is that the principal symbol of the Pauli-Fierz 
Hamiltonian, defined in (fT^ . is diagonal with respect to the Fock projectors 
Qmi while for the Nelson Hamiltonian one needs a dressing transformation 
already at the leading order. This makes the analysis of the Pauli-Fierz case 
somewhat less technical. 

The effective dynamics for M = (dressed vacuum) was calculated by 
Spohn ( |Sp| , section 20.2) in the case when the photon has a small mass, 
> 0, which introduces a gap in the principal symbol of the Hamiltonian. 
He however states that these effective dynamics are identical with the ones 
calculated for the case c — > oo, while we have already remarked that the 
spin dependent term cannot be present when e — 0. In the case mph > 0, 
moreover, the transitions between the almost invariant subspaces become 
smaller than any power of e and therefore it is not known how to get an 
explicit expression for them. 

The Pauli-Fierz Hamiltonian has also been extensively studied to get 
informations about its spectral and scattering structure. Not pretending 
to be exhaustive, we refer the reader interested to these aspects to |BFS] . 
|DeGej . |FGSj . |GLLj and references therein. 

In section [2] we complete the description of the model and discuss the ap- 
proximation of the original dynamics through infrared cutoff ones. In section 
[3] the construction of the dressing operator ^ is discussed, and applied in 
section H] to the study of the dressed Hamiltonian. The main results on the 
effective dynamics and the radiated piece are contained in section O Finally, 
appendix [Xl contains a sketch of the proof of theorem [H for the case c — ^ oo. 

2 Preliminary facts 

In this section we elaborate on the definition of the Pauli-Fierz model and 
discuss some preliminary facts like the self-adjointness of the Hamiltonian 
and the approximation of the original dynamics through infrared cutoff ones. 

2.1 Fock space and field operator 

(The proofs of the statements we claim can be found in ([ReSi^j, section 
X.7)). 
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We denote by the subspace of the Fock space, defined in ([7]), for 
which = for all but finitely many M. Given / e L'^{^^ x Z2), one 

defines on the annihilation operator by 



(a(/)^)(^)(A;i,Ai;...;fcM,AM) :=v/mTTV I 

A=i 



dkf(k,xy- 

^ ^ (35) 



■¥'''+'\k,X;h,Xi...,kM,XM). 

The adjoint of a(/) is called the creation operator, and its domain contains 
<^{in- On this subspace they satisfy the canonical commutation relations 

[a{f),a{9T] = (/,^)l2{r3xZ2), , . 

[a{f),a{g)] = 0, HfY ,a{gr] = . ^ ' 

Since the commutator between a{f) and a{f)* is bounded, it follows that 
a(/) can be extended to a closed operator on the same domain of a(/)*. 
On this domain one defines the Segal field operator 

m ■■= ^w) + aim (37) 

which is essentially self-adjoint on ^fin- Moreover, ^fin is a set of analytic 
vectors for $(/). From the canonical commutation relations it follows that 

[^f),^g)] = iQ{f,g)L^u^^^,)■ m 

Given a self-adjoint multiplication operator by the function u on the 
domain D{uj) C L^(R^), we define 

^^,fi„ := C{Qp, a(/i)* ■ • ■ aifMT^F : M G N, /, G D{u)^C^ j = 1, . . . , M}, 

(39) 

where C means "finite linear combinations of". 

On ^ujfin we define the second quantization of u, dr(u;), by 



M 

idT{u)^Y^'\k,, Ai; . . . ; ^m, Am) := a;(A;,)^(^^)(A;i, Ai; . . . ; ^m, \m) 

dr(tu)nF := 0, 

which is essentially self-adjoint. In particular, the free field Hamiltonian Hf 

acts as 



A/ 

(i7f^)W(A;i, Ai; . . . ; ^m, Am) = 5^|fc,|vI/W(A;i, Ai; . . . ; kM, Xm) 
Hi^F = 0, 
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and is self-adjoint on its maximal domain. 

From the previous definitions, given / G D{uj) ^ C^, one gets the com- 
mutation properties 

[dlM, aim = aiuJ fT, [dlM, a(/)] = -a{uf), 
[dr(a;),i$(/)] = $(ia;/). ^ ^ 

2.2 The Pauli-Fierz model 

Using the Segal field operator one can write the quantized vector potential 
and the magnetic field appearing in ([H]) as 

A^{x) = <l>{v,), (41) 
v,{k, A) := f{k, X)e-''-^, f{k, A) := ^j^^0ik), (42) 
B^{x) = V^. X A^{x) = -<l>{ik X v.^), (43) 

where ex{k), A = 1,2, are, for simplicity, real photon polarization vectors 
satisfying 

ex{k) ■ e^{k) = 5x^, k ■ e^ik) = . (44) 
The smeared Coulomb potential is given by 

V^^coui(x) = ^ /" dA;e^^-(-^--')^^ . (45) 

Analogous expressions hold for the infrared cutoff Hamiltonian H'^''^, where 
the form factor ip is replaced by (pa- 

To separate more clearly the terms of different order in the Hamiltonian 
H^, equation (ITTl) . it is useful to write it as 



H^ = Y,e'hl (46) 



i=0 



where 



^ 1 

N 



hl--=-Y.±f'rHv^,l (47) 

i=i ^ 
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Each of the hi is of order 0{1) when apphed to functions of bounded kinetic 
energy. The coefficients for H^''^ will be denoted by /i^'*^. 

As proved by Hiroshima [Hi] using functional integral techniques the 
Hamiltonian (and analogously H'^) is self-adjoint on Jifo for every value 
of the masses, charges and number of particles. Since however we study the 
limit £ ^ (respectively c — oo) it is enough for our purposes to show this 
using Kato theorem, like, e. g., in [BFSj . 

Even though the proof is well known, we repeat it because we need to 
show that the graph norms which appear are equivalent uniformly in e and 
a. Moreover, the estimates which appear in the proof will be useful in propo- 
sitions [3] and lemma [H 

Given f e L'^ {R^ x Z2) , we define 

ll/L := (ll/|A:r^/li.(M3xZ.) + ll/lli^(M3xz.))'/' ■ (49) 
One has then the basic estimate 
Proposition 1. 

■ ■ ■a"(/n)(i^f + < all/lL ■ ■ • WfnL , (50) 

where a^{f) can be a{f) or a*{f). 

Proposition 2. Both Hamiltonians and H'^''^ are self-adjoint on Jifo. 
Moreover the graph norms they define are equivalent to the one defined by 
Hq uniformly in e and a. The same holds for the graph norm defined by 
{H'f^ and {H^'^^y/^ 

Proof. (We give the proof for if^''^, the one for H'^ is the same) 

The regularized Coulomb potential is a bounded function, therefore for 

it the statement is trivial. 

We choose a vector in a core of made up of smooth functions with 

compact support both in x and k. 
For the term of order e we get then 



N 

e 



j=l J 



< 



N I I 3 N , ^ 



^ — ' rrij ^ — ^ ^ ^ — ^ rrij ^ — ^ ^ 

j=l ■' a=l j=l OL=l 



< 



^3 

N I I 

e 



^-^ rrii 
j=i J 
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so for e sufficiently small this term is Kato small with respect to the free 
Hamiltonian, with a constant uniformly bounded in e and a. 

An analogous estimate holds for the term with the magnetic field. For 
the remaining one we have 



-2 ^3 



N ^2 



2m- ^ " ' 



N ^2 



< 



which completes the proof. 
Proposition 3. If (y{e) < then 



c{M,.^) ^ C\t\a 



1/2 



□ 



(51) 



Proof. From the previous proposition we know that both Hamiltonians are 
self-adjoint on J^, so, given ^ G Mq, we can apply Duhamel formula to get 



Putting := e ^^^'/^^i ^ the difference of the two Hamiltonians is 



AT 



N 

E 



j=i J 



(52) 



Using the explicit expression ( H5l) . the term with the Coulomb potential 
gives 



ipa COul 



x)\ < 



1 ^ 



r dk 


1 - |2 1 " |2 


1 \k\^ 





127r 



1 ^ 

W2 X] 
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=^ IlK^coui — V^^cou\\\c{.jr) = ^i^r) ■ (53) 

For the term of order e, proceeding as in the proof of proposition [2] we get 
that 



AT 



j=l J 



N 



< ^||l(o,a)(/^)t^x,,<7LI|^slUo 



j=l 3 



From the same proposition it follows that the graph norm associated 
to Hq and the one associated to H'^ are equivalent uniformly in e and a, 
therefore 



N 



N 



(54) 



The same reasoning holds for the term containing the spin, which has 
however a. \k\ more, which gives in the end 



N 



E • ^^-'^ ^ l(0,.)(fc)t^x,)^. 

i=i ^ 



0{e'a'/')mU ■ (55) 



Concerning the last term we have 



N o 



N o 



2m i 

=1 



2m 
j=i 

N ^2 



^'I]^[*K'-)'^(l(0,-)(^)^-.)+'^(l{0,<x)(^)^x,)'^'(^^x„.) + 



2m 

+ <l>(l(o,a)(A;)i;xJ^]^s . 
Using again the basic estimate in proposition [1] we get, for example, 

\\^{Vxj,(j)^{'^{0,a){k)Vxj)'^ s\\ < C||^2;^,o-|U||l(0,o-)(fc)t^Xj|| 

by the same reasoning we used for the terms of order e. 



(56) 
□ 
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Lemma 1. Given a function x G C^(M) and assuming a{e) < e^, then 

\\xm-xiH^n\\ciM^)<Csa'/' (57) 



Proof. Using the Hellfer-Sjostrand formula (see, e. g., |DiSj| chapter 8), given 
a self-adjoint operator A, we can write 



X{A) = - [ dxdydx%z){A-z)-\ 



z := X + iy, (58) 



where x° ^ C'^(C) is an almost analytic extension of which satisfies the 
properties 

ViVeN 3Dff : iBx^l < Dffl'^zf, 
X|b — X 

(For the explicit construction of such a see |DiSj| ). 
Applied to our case (13^ yields 

Xm - xiH^n = - [ dxdy dx^iz) [{H' - z)-' - {H^^^ - z)-^] . 

Since both Hamiltonians are self-adjoint on J^q we have 

{H' - z)-^ - (i/"'" - z)-^ = {H'^" - z)-\H'''' - H'){H' - z)-^ , (59) 
and hence 

\\x{E^)-x{ii^n\WM^)< 

< - [ dxrfyiax"(^)iii(i^^''^-;^)-i£(,.^')ii(^''''-^')(^'-^)"'iU(.^)- 

In addition we have that 

UH^'''-z)-^ci..r)<j^. (60) 
This follows because is dense in the domain of ff'^=o,f^=o^ ^nd for every 
H^,<^^ H'=^'^=^^ as (e, a) (0, 0) . 



According to theorem VIII. 25 |ReSii| , this implies that 
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therefore I'^zlW {H'^''^ — z)^^'$\\ is bounded for every \E' and the uniform bound- 
edness principle gives f lUUj) . 

For the second norm we find that for z G supp 

\\{H''- - H%H' - z)-'\\ci.^) < 

< II (i^^''^ - H'')\\c(M,J^') ■ II (-^^ - 2)"l£(.^V^6) 

The right-hand side was aheady estimated in the previous proposition 
(see equation ( 152|) and the following calculations), the only difference being 
that here we have to multiply the final result by £. □ 

3 Construction of the dressing operator 

In this section we construct the unitary dressing operator ^g^^. We apply to 
the Pauli-Fierz Hamiltonian the general procedure explained in some detail 
in |TeTej and, for the case with spectral gap, in All the calculations 

expounded in section [3?T] are formal, and they serve as a guide for the rigorous 
definition of ^e,a given in section I3.2[ 

3.1 The formal procedure 

The main idea is to build an approximate projector, 7r*^^\ which satisfies 
formally 

(^(l))2 _ ^(1) ^ ^^^2)^ [7r«, if^] = 

Integrating over time the second equation one gets in a loose sense that 
[g-iiHVe^^(i)] = 0{e\t\). 

The projector 7r(i) is found using an iterative procedure, which assumes 
that one can expand it in powers of e, 

vr*^^^ =71-0 + en^, 

where the coefficient ttq is a known input and must commute with the co- 
efficient of order zero in the expansion of the Hamiltonian H'^, see equation 
(HUj) . As it turns out, the procedure does not work directly for H'^, but only 
for the infrared cutoff Hamiltonian H'^''^. 
An obvious choice for is 

^0 — Qm, 



18 



which satisfies [11^'^ ,Qm] = 0. 

Proceeding now in the same way as described in [TeTej we get a formal 
expression for the first order almost projection given by 



7I"m '■= Qm + 

N 



M 



M 



^-^ rrij 



ifz,-,a(A, k) 



\k\ 



For brevity, we put from now on 

fwx^^a{\, k) 



so 



It is clear from equation fl6ip that 

QmtTi^Qm = (1 - <5M)7rj^(l - Qm) = 0, 

fc^^j is also almost invariant for the total dynamics, in the sense that 
[7r£), = [Qm, K'''] + e[7rf, h^^ + sIQm, + e'[n^ , hr] + 

To justify this claim we note that 

N N 

[^f , Kn + [Qm, hr] = iY.^[pr [Qm, K-] - E ^p,. 



(61) 



AT 



rrii 
■',1=1 J 



[Qm, = i 5Z — [Pj ■ [Qm, y] + 



AT 



N 



i — [Pj ■ [Qm, ^j,ai V^^ coul] +'^Y — [Pr [Qm, Hi] + 



N 



TV 



I I L r, I I L It 
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To analyze in a simple way the restriction of the dynamics to the subspace 
defined by tt^' one builds an almost unitary U^^\ which maps the almost 
projections to a reference projection up to terms of order 0{e'^). Using the 
formal expression we get for we will define in next section a true unitary 
operator which will allow us to construct a rigorous version of the almost 
projections tt^^ . 

A natural choice for the reference projections, hnked to the physics of the 
system, is to choose them equal to the QmS. We assume then that also U^^^ 
can be expanded in powers of e, 

U^^^ —1 + eUi, 

with the condition Ui + — 0{e). This ensures that 

To determine Ui we impose that intertwines the almost invariant 
projections with the reference projections Qm up to terms of order 0{e^): 

The left-hand side gives 

= Qm + Qm] + vrf ) + = QM + e [[U^, Qm]+ 

N 



3 

SO we can choose 



1 \ 

,=1 "^3 \ J 



N 

e 



^^'^ = l + i^E^^^>-^^- ■ (62) 



j=i -I 

3.2 Rigorous definition 

To get a well-defined unitary operator from the formal expression for U^^"^ we 
first cutoff the number of photons in the field operator ^j,a, replacing it by 

:= g<L$,-,.Q<L, (63) 



where L is fixed, but otherwise arbitrary. 
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We introduce then a cutoff in the total energy, to cope with the unbound- 
edness of the momentum of the electrons pj . This reflects the fact that the 
adiabatic approximation holds uniformly only on states where the kinetic 
energy of the slow particles in uniformly bounded. 

More precisely, given a function x £ C'^(M), we define 

N N 

3=1 ^ i=i ^ 

= 1 + ExiH^nut + ^(1 - x{H^n)utx{H^n, 

(64) 

where we have defined 



N 




Note that it follows from the Coulomb gauge condition that U]^ * — —U\^^^ 

so 

U'il; = 1 - exiH^nut - e{l - x{H^n)Utx{H'n ■ (66) 

Lemma 2. 1. For each x e R^, ^]^^^{x) e and $};^(a;)* = 

Moreover, 

^l,:R'^C{^), x^^l^ix), eCb°^(M^£(^)), (67) 
and, for a small enough, 

W^lJciJ^) < CVLTTVlog(a-i). (68) 

Given a with \a\ > 0, it holds instead 

d:<^la = d:<^lo + 0(al"l VL+T)z:(^), (69) 

where 

is a well-defined bounded operator on M' . 

2. The statements of point 1 (except for the self-adjointness o/$^^(a;)j 
remain true if ^ is replaced by D{Hi). 

Corollary 1. The fibered operators 5^$^,^ belong to £(^)n£(^) Ma e N^. 
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Proof. The proof of both statements follows from the facts that 

\\Q<LHgx{-))Q<L\\c{.-r) < 2 1/ Viv + 1 sup 1 1 ( ■ ) I U2 (K3 X Z2 ) 

and that 

d:v^{k,X)\k\-' = {-ip\k\~^^^-^\^\ex{k)0{k)e-''-^ . 



□ 



Lemma 3. The operator U^^^ is closable and its closure, which we denote by 
the same symbol, belongs to C{J^) fl C{J^o)- Moreover 

mi'iwcij,) < c(i + £Viog(a-i)), (71) 

where = M' or The same holds for U^^^* . 
Proof. The operator 

is defined on D{j>j) and, since x{.H'^''')^\(t is a bounded operator, we have 

which is clearly bounded. This shows that -pj is closable and its 

closure belongs to C{Jif). 

The same reasoning can be applied to the operator 

which shows that U^l^ is also in C{J^q). 

The estimate on the norm follows now from the estimate on the norm of 
"^jcr given in lemma [21 □ 

Theorem 2. Assume that a = cx^e) satisfies the condition 



£v/log(a(£)-i) 0, £^0+, (72) 

then the operator 

^ ■■= (73) 
is well-defined and unitary, for e small enough. 

Both and belong to C{J^) fl C{J^), with the property that 

W'^WdM), \\^*\\c(M) <C, (74) 

where C is independent of e and a. 

Moreover we can expand them in powers of e and the corresponding series 
converges both in C{J^) and C{J^q). 
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Proof. It follows from equations flM|) and flUBj) that, defining 

:= xiH^nUt + (1 - (75) 

we have then 

f/i;^ = 1 + £T^(,), T; = -T,, \\T4c(.x) < Cv/log(o-i), 

where ^ = ^ or J^q. 

From this expression we get immediately that 

SO, choosing e small enough, we have that e^HT^^^^ ||£(^) < 1, therefore the 
square root is well-defined, and can be expressed through a convergent power 
series in C{J^): 

From standard calculations it follows in the end that is unitary on 

4 The dressed Hamiltonian 

We define the dressed Hamiltonian as the unitary transform of H^'" , 

H^f^^ := ^H''"^* . (77) 

Since is a bijection on -ff^res self-adjoint on and using the 
expansion of on C{J^q), we can expand -ff^res ^^i.'^Oi ^)- 

Theorem 3. Assume that a = cr{e) satisfies conditions (1251) . The expansion 
of the dressed Hamiltonian up to the second order is then given by 

where 

^ 1 

/^r = ^0 = E + ^^--1 + 
j=i 

hf^"^^ is given in equation (IHTI) and h'^'^^ in equation flH^ . 
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Proof. Applying equation flTHj) . we get that 

^ = 1 + eT + — + 0{e^T^'Q, 

therefore 

i/^;^^ = (1 + + e^T''l'l)H^'\\ - eT + e^T^ /2) + 

= hr +e{hr + x{H^n[ut^ m + "(i - x) ■ • • 



=: + ehf^' + ^'/^f ' + (^(^'(logo^"')'/')^^„,^) 

where "(1 — x) ' ' "X" indicates that for every term containing x ■ ■ ■ we have 
to add a corresponding term containing (1 — x) ■ ■ ■ X; as in equation (|75il . 
Using equation fl53l) to ehminate the a, we get immediately 

ht'' = h ■ (79) 
The commutator in the term of order e gives: 

N 

[f^t, ^ol = -Q<L/^rQ<L + e — *t ■ V-.^^^ COU1 + 

(80) 



therefore, taking into account equation 

^res ^ (1 _ ^(^H'^<^))Q^^hlQ<^{l - Xl^'"'^)) + Q<LhlQ>L + g>L/i!g<L + 
+ Q>L/ifQ>L . 



^The expansion of till the second order coincides with that of e'^'^ , however e^"^ is 
not the correct dressing transformation to every order. Following the formal procedure 
sketched in section [01 one can construct a second order expression for , which depends 
however also in general on the second order coefficient of the Hamiltonian, h'^'" and has 
not a simple exponential form. 
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We will show below that this term vanishes in the effective dynamics. 
Concerning the terms of second order we have 

N 

Pt, hin = -i E ■ *(-^.^) ■ pi] = 

„ (82) 

3,1=1 ■> 

where 7?.l-i is a term which vanishes on the range of Qj when j < L — 1. 

Using equation ( l38i) to calculate the commutator of the two field operator 
we get in the end 

- ' E ^r'- ^"^'^^ - E ;^«(ff • 

j,Z=l -' j,'=l 

= i y /■ rffc M^e^fc-(-.-Op . . (1 - ® = 

= i y / M^e^'^ -^p . . (1 _ ® K)p,e-''=-^' 

The remaining term of second order gives, 

[T^xiH^nPt^hr]] ^ -[T,x{H^nQ<LhrQ<L]+o{s) ^ 

= -xPt' x]Q<LhrQ<L - x'Q<L[Ut, hr]Q<L - Xix, Q<i.hrQ<M, 



^3) 



Putting the calculations above together, and using equations ( !55l) and 
561) . we get in the end 



N N 

ht' = k + Yl — "^t ■ v., V^coui + Yl ^^V., • P,) ■ Pi+ 
+ XPI.. /^I'l - ^x'Q<L[t/t' ^i"]Q<L + (1 - X)Pl.^ ^^^^ 

- \x\x.Q<ihrQ<^Pl. - \xPl„.x\Q<ihrQ<i^. 

where the commutator is given in equations (!H2|) and (!H^ . We 

will show below that only few of the terms written above contribute to the 
effective dynamics, giving the correction to the mass of the electrons and the 
Darwin term. □ 
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5 The effective dynamics 

We quote without proof a number of lemmas, which, with minor modifica- 
tions, are identical to the ones proved in [TeTe] . 

Lemma 4. (see corollary 4 ITeT^ ) 

Given a function x ^ C'cTl-'^) '^^'^ a a > 0, we have 



where ^ G C^(M), = X one? 



where 



:= 2c^ + i?oo + minjc;^. A} 



So 



sup{\k\ : G supp x} 
sup |Kpcoui(a;)| , 



^5) 



and we can choose sup{|fc| : k G supp ^} arbitrarily close to Cg. 
74n analogous statement holds for the creation operator. 

Lemma 5. Assume that a satisfies conditions (1251) . then 

1. Given a function x G C^(M), we have 

xiH^n - Xik) = ^^x' (86) 
where 71^ G C{J^,J^o), = C'(l) and 

n'^QMXiH) = (Qm+1 + QM^imH)n'^QMxChl) + O{e)cij^,jroh (87) 

where C, has the properties described in lemma^ 

2. Moreover, we have that 

mZs) - xik) = c?(^)£M), (88) 

and that 

Qmx{HZs) = QMximWL) + (^(^'v/M^)£M), (89) 
where x is any C^(M.) function such that xx = X ^^^'^ XX = Xj M < L — 1. 
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Proof. We give the proof just for point 1. The proof of point 2 is analogous 
and can be found also in ( [TeTej . lemma 7). 

Proceeding as in lemma [T], we apply Hellfer-Sjostrand formula and we get 

xiH^'^) - Xik) = - [ dxdydx^iz)[iH'^- - z)-' - CK- z)-']. 

Since both Hamiltonians are self-adjoint on the same domain we get, 
iterating equation fl3I?|l . 

{H'''^ - z)-' - Chi - z)-' = -sCK - z)-^hr(hl - z)-^ + 0{e''\^z\-'')ci.^,^,y 

The statement now follows from the explicit expression of h^{^ using 
lemma HI □ 

Theorem 4. (Zero order approximation to the time evolution) 

||(e-i^r4 -e-^o!)Q^,;^(i7^;-)||^(^^) = 0(v/MTl|tkVlog(a(£)-i)), (90) 
||Q,,(e-i^r4 -e-'^o|);^(i7^;^j||^(^.) = 0{^MTl\t\e^\og{a{e)-^)), (91) 
for every x G C^(M) such that XX = X- 

Proof. Using lemma [5] we replace x(-f^dres) with x(^o)' since the difference, 
being of order 0{e), is smaller than the error we want to prove. 

Both Hamiltonians are self-adjoint on J^, therefore, applying Duhamel 
formula, we get: 

(e-^*^-^=/^-e-'*^o/^)QMX(/^o) = 

= -- / c?se'(^-*)^^^^.A(i7^;- _/,g)e--^§/-Q^ 
^ Jo 

= _i /■ rfse'(^-*)^^^^e/^/ii,;,e-^^'^°/^gMX(/i^) + 0{e^/WTl^\og{a-^))ci^) 
Jo 

= _i /■ rfse'(^-*)^d;:e^/,fesg^~^^.)g-is/^o/. + 0{e^W^^\og{a-^)) c^^^y 
Jo 

Putting in the previous equation the expression of h'^'^^, equation (IHTl) . and 
observing that we can replace, again by lemma [5l x(if^'°") with x(^o)? S^^ 
that the right-hand side is of order 0{\/M + l\t\e^J\og{a~^)). 

For the second estimate we apply again Duhamel formula, inverting the 
position of the unitaries, 

QM(e-*^-/^-e-'*'^o/^)^(^^^-) = 
^ Jo 

= _i rcise'(^-*)'^o/^gM/^^^x(i^d;:s)e-'^''--/^ + 0{E^MTl^\og{a-')) , 

^0 
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so, replacing xiH^ms) xi.H^''^) by x(^o)' claim is proved. □ 

Corollary 2. The dressed projectors Plj are almost invariant with respect 
to the original dynamics, i. e., 

\\[e-''''^,P!,]xm\\ci.^^) = 0{^MTl\t\e^/\og{a{e)-^)) . 
Proof. From the previous theorem it follows that 

\\[e-''''^-^',QM]x{HZs)\\ci.^) = 0{VMTl\t\e^/\og{a{ey^)) . 
From the definition of -f^^^gg, equation fl77|) . we deduce therefore that 



but 



[e-'''''"^,Plr]x{H'n\\ci.^^) = 0{VMTl\t\e^\og{aie)-^)) , 

+ 0(v/MTT|t|£Vlog(a(£)-i))^(^), 

where we have used equations (ISTl) and ( |571) . the fact that xi.H'^) ^ Mq) 

with norm uniformly bounded in e and equation fl74l) . □ 



Lemma 6. The truncated dressed Hamiltonian 

HZ:=hl + eht^^^ + e'M^^^ (92) 
is self-adjoint on for e small enough. 

Proof. The proof follows from a symmetric version of the Kato theorem 
( IReSial , theorem X. 13) . □ 

Theorem 5. (First order approximation to the time evolution) Given a 
function x G 



e-''"^-'^QMX{HZ.) = 

= e-''</^QMX{HlfJ - t ds e'(^-*)^o/^/i2,ODe-'^^o/^QMX(i^d';:s) (93) 

Jo 

+ 0{e'/^\t\)ci^^{l - 5mo) + 0[e^^J\og{a-^){\t\ + \t\^))ci.^^ , 
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where 6 mo = 1 when M = 0, otherwise, [Hj^ ,Qm] = WM, 

j=l J 



AT 

_2 " 



^^^m^m^iKS 2 1 A; I 
The off-diagonal Hamiltonian is defined by 



N 

j=i ^ 



h2,OB ■■= —^ho ■ Va^.K^coul • (95) 



Remark 4. The spin term, as mentioned in the introduction, does not appear 
in the effective dynamics, even though it is apparently of order 0[e^). This 
is due to the fact that it is off diagonal with respect to the decomposition 
of the Hilbert space associated to the Qms and that the coupling function 
in the magnetic field, equation ( H3|) . goes to zero like when \k\ — >■ O"*". 

This implies that the term is actually smaller than 0[e^\ as explained in 
the proof. 

Proof. We split the proof into three parts. In the first one, we show that 
equation (!93|) is true with a diagonal Hamiltonian H^^'' given by 

N N 2 



i=i ^ i=i ^ 



(96) 



i,j=i ^ 

and an off-diagonal one /i2,oD defined by 



^ e • ^ e ■ 



j = l J j,l = l J 

N N 2 

In the second part we prove that if one neglects the term 

N g2 

^'E^«K)*«K) 

j=i J 
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(97) 



in , one gets an error of order in the time evolution. Note that 

this term is exactly zero if the initial state for the field is the Fock vacuum. 
In the third part, we prove analogously that we can replace /i2,oD with 

^2,OD- 

More specifically, the terms which we neglect in Hj^ and /?.2.od give rise 
to higher order contributions to the time evolution, although their norm in 
C{J^o, is not small. 

This is caused by the fact that they are strongly oscillating in \k\, so 
that their behavior is determined by the value of the density of states in a 
neighborhood of A; = 0. For all these terms, the density however vanishes 
for k = 0, uniformly in a, and this implies that they are of lower order with 
respect to the leading pieces whose density is constant (for the terms in H^^^) 
or diverges logarithmically in a (for the terms in /i2,oD)- We elaborate on 
this last observation in a corollary to this theorem. 

We start showing that we can, up to the desired error, replace -ff^res 
-^dres- -^y lemma [6l H^^l^ is self-adjoint on J^q like -ff^^es' therefore we can 
apply the Duhamel formula and use theorem [3] to get 

e-i*^d;:e/^ _ e-'*^^l/^ = -1 /*rfse^(^-*)^^^^=/^(ff^;:3-i/(l)e-^^^^l/^ 

^ Jo 

Moreover, using lemma El we can replace QMx{H^fes) by QMX{ho)x{H^fes)- 
Since the diagonal Hamiltonian if^ is also self-adjoint on for e suf- 
ficiently small (the proof can be given along the same lines of lemma E]), we 
apply again Duhamel formula, 

^ Jo 

= _i / tisei(^-*)^dl/^/ifeSg-is/f(^)/.g^^~(^.)=(^.,.j 

Jo 
Jo 

To analyze the first term, we remark that, proceeding as in lemma El one 
can prove that 

so 

[e-<'/^x(/^g)] = O(e^v/M0)£(.ir), 
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therefore, 

-i / ds e'(^-*)^^?^e/^/ir^e-'^</^gMx(/^g)^(i^d;:j = 

Jo 

+ 0(£2|t|v/log(a-i))^(^). 
From equation flHTj) it follows 

/^f^QAfx(/^^o) = (1 - x(i^^"^))Q<L/^IQ<L(i - x(i^^''^))QAfx(/^g) = 

using lemma [5] twice and lemma HI 

Concerning the second one, applying once again the Duhamel formula, 
we have 

^0 

+ 0(£^|t|Vlog(^-^))£(.i^), 

(2) 

so we have to look at e^'-*~*''^drosA(/idres_^^ .^-jq^^^^^e-j^ Following a procedure 
already employed several times, we first observe that, in the expression for 
hf^^^, equation (18^ . we can replace, making an error of order 0{e), x{H^''^) 
with xCK)- 

Applying lemma IH we see that the last three terms in (1841) vanish, and 
that the terms containing [Ui^, h\''^] combine to give the Darwin term and the 
expression of the effective mass. What remains is exactly /i2,OD- Finally, we 
apply again Duhamel formula to approximate the time evolution generated 
by H^iL getting 

= e'^'^-'^^<^/%oBQMxCK) + 0(£|t|Vlog(a-i)) . 

~ (2) (2) 

We proceed now to show that we can replace if^ with H^^ , up to an 
error of order 0{£'^^'^\t\)c{jf)- 
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Applying repeatedly Duliamel formula, and putting \1' := QpixiH^ZsJ^^o, 
we get 



N 2 ft 

+ 0{e^\t\'')ci^). 



,=1 -'^Jo 



To streamline the presentation, we assume that M = 1, the calculations for 
M > 1 are basically the same, but more cumbersome. 
The integral gives therefore 

e-^*'^o/(A;i, Ai) V / dkf{k,Xy [ c/se't(l*^il-l^l)e't''pe'"^-(^-^^)e-'t''p^(fc. A) 

ilR^ Jo 

= e-i%(A:i,AOX^ / . ^ -/it ^ ill V-i + " 

^ V 1 + i(|fci| - 1 

■ rdse't(l^'il-l'=l)e'f''pe^"^(^-^^)e-^t%^(A;,A). 
io 



(9^ 



Integrating by parts we get 

i{\ki\ - \k\)e~'^ [ cise'"(l'=il-l'=l)/V*''pe^^^-(^-'=i)e-^'''p^(fc,A) 
Jo 

_ ^it{\ki\-\k\)/E^ith^^ix.j-{k-ki)^-ith^-qj _ gi3::j-{fc-fci)^_|_ 



- [ dse''(l'=il"l^l)/V'''p[/i^,e'^^-(^-*^^)]e-'''^p^, 



where the commutator is of order 0{6) when applied to functions of bounded 
kinetic energy, so that the right-hand side is uniformly bounded in e. 

We have now to put this expression back in ( l98l) and estimate the single 
terms. We show how to do this for the first one, the others being entirely 
analogous. We ignore the unitary on the left, which does not change the 
norm, so we have to consider 



f{k,xy 

\k\)e-^ 
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Using twice the Cauchy-Schwarz inequality we get 

2 



|2 



< t 



ai,...,(Tff ,Xi = l 
2 

E 

(Tl,...,CTJV,Al=l 

2 . .t 



2 

da; /" rffci|/(fci,Ai)|2^ /" dk- 

A=l 



l/(fc,A)p 



+ {\k^\ - \k\ye-^ 



A=l 



dk ds 



\ \ 1 



Ai,A=l 



| /(fc,A)/(fci,Ai)p 
+ (|A;i|-|A;|)%- 



< 



dki 



A/e 



dk 



kik 



1 + (A;i - A;)2 



.if'. 



We proceed now to examine the last three terms in ( 1971) to show that 
they can be neglected. 

The second and the third term are a sum of terms of the form 

^{giX,k)e-''-^^)f, 

where the function g ^ \k\ —>■ O"*", with a = —1/2 or +1/2, and T is a 
Pauli matrix, or the product of two momentum operators. 

For a term of this form we get then, putting again \E' := (5AfX(-f^drcs)^0! 

ie f ds e'(^-*)^oA$(^(A,A;)e-i'=-^^)fe-'^^o/-g^,;^(if^-;-Jvl/o = 



= i^ / ds e^(^-*)^o/"[a(^e-^^-^0 + a{ge-''^''^y]fe-''^y'^> . 
v2 Jo 

Expression of this type have already been estimated in ( |TeTej . theorem 4). 
For convenience of the reader, we give the proof for the annihilation part, 
referring to [TeTej for the creation part, which is entirely analogous. 
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The annihilation part gives 

2 

i^^e 



v2 jQ f— 

(99) 



A=l 

■ ^(x, a; fc, A, fci, Ai, . . . , /cm, Am 

2 



• /" rfse-''l^l/V'''p/V'=-^^fe-^^'^p/"^(x,a;A;,A,...) . 
Jo 

Integrating by parts we get 

- i\k\e-^ [ ds e~''m/eQishl/e^.k-x,f^-\shl/e^(^^^ ^. f.^ X, . . .) = 

Jo 

-L f rfse-^'l'=l/V"'^p/^[/i;,e^'=-^^f]e-*"'^p/^^, 
^ Jo 

where 

[hl,e''-^f] = -L(2ee''-''^k-p,+e^\k\'e''-^^)f + e'''-^^[^^^^ . 

The commutator on the right-hand side is zero if T is a Pauh matrix, or it is 
a term of order e times p, if T is the product of two momentum operators. 
Therefore it has the same form as the first part of the right-hand side and 
can be treated in the same way. 

We have now to put the result of the integration by parts back in equation 
flU^ and estimate what comes out. We show how to do this for the first term, 
all the other ones can be treated in the same way. 

Ignoring the unitary on the left and the constants, we consider then 

eY f dk Trfse-'^l'^'/V^'^p/V'^-^^fe-'^'^p^^ 
fr'^jRS l-i\k\e-^Jo 

Using the Cauchy-Schwarz inequality we get 

■\t\Y f dk [ ds|e^'''p/V'=-^^Te-''^p/"^|^ 
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so 



A=l 

The left integral gives 

12 /-A 



Jr^ -L "t I'f'l ^ Jo 



when a = 1/2, or —1/2. 

The same analysis can be carried out for the remaining term 



The proof is identical to the one given in ( |TeTej . theorem 4) for a similar 
term appearing in the case of Nelson model, and depends as above on the 
fact that v^^{X,k) ~ \k\-^/^, \k\ 0+. □ 

Corollary 3. At the leading order, the radiated piece (i. e. the piece of 
the wave function which makes a transition between the almost invariant 
subspaces) for a system starting in the Fock vacuum, "^olx) = iP{x)Qf, 
ip{x) G is given by 

(100) 

where xf is the solution to the classical equations of motion 

mjxf{s;x,p) = -Vx,V^couiix''\s;x,p)), ^^^^^ 
xf{0; x,p) = Xj, xf{0; x,p) = PjmJ^ j = 1,...,N. 

This coincides with the leading order of the radiated piece corresponding 
to the original Hamiltonian H"^ , for a system starting in the dressed vacuum. 
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Proof. Applying equation fl93|) for the case M = we get at the leading order 
Q^e-'*^d;:sA^(a;)^]p = _i£ f ds e'(^-*)'^o/"/i2,oDe-"^o/"V^(x)(]F = 

v2 ^"^j io 

= -^^^^(^) ■ £ ^ / e'(^-)IWv(^-)%A(e--. - 1) 



+ 0(£'|t|)£(,ir)||V^|U2(R3n), 



where we have used Egorov's theorem to approximate e'^'^p/'^V^jV^couie '^'^p/^ 
(see, e. g., [Ro]). 

To end the proof of the first statement we have to show that the norm of 
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the first term is small. For the jth term in the sum we get 



spin 



dse^'l'^l/V'p/"(e''=-^^ - !)■ 



■Vxjy^coui(x)e' 



X] 



J spin 



/ c?x dk 



\k\' 

2 



< e^\t\^ I dk 



2 ft 



ds 



= Ce% lo; 



A/e ^ 

c/|A;|— / ds 

C7(e)/€ 



2 



A 



ds 



^i^)/ Jo 
For the second statement we have 



1 - Po^)e-i*^^Po^x(i^^)^ = (1 - n')(e- 



e ^ 



+ (1 - Po')e-^*^ P'^\x{m - xiU'^^W + (1 - ^o)e-^*^ ' Plx{ii'n'^ 

= Qo^e-'*^^^ce^^(x)^]F + 0(e'loga(£)-ivI/||^^.) , 
with 

Remark 5. Denoting by 



□ 



xf{s]x,p) 



the operator acting on which appears in (llOOp . the norm squared of the 
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leading part of the radiated piece is 



e 
~2 



—in A=l,2 J J II 



^crfek-i I'^l 



'(T{e)e 

where a > can be chosen arbitrarily small. 

The symbol of Tt{k) is an e-independent function, which for = is 
different from the null function 

so Tt{0) is different from the zero operator. We expect therefore that for a 
generic state ip 

inf ||T,(|A;|)^||,^^^c3 > . 

0<|fc|<a 

In this case one gets as lower bound for the norm of the radiated piece 
inf ||T,(|A:|)^|U^^c3^f T ^1^1^) ^ = 0(eVlog(ea(e)-i)), 

which is almost of the same order as the upper bound. 

Remark 6. The radiated energy, defined in equation ( l33l) . can be written at 
the leading order as 

£^rad(t) = (Qo^e-'*^d;:e/-^(a;)f]p,i7fQ^^e-"^d;:s/-^(a;)^]p), 
where ip is defined in (I3T1) . Using the expression for the radiated piece we 
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get 



t/|A;||^.(,)(|A:|)n|T,(|A:|/e^(x)||^^^c3 = 
d\k\\0{\k\)\^ f ds f ds'e^^'-''^\^\- 



2 



spin A=l,2 ■ 

3' 
4 
3' 



-7r£: 



^0 



where we have used the product formula for pseudodifferential operators (see, 
e. g., |Roj ) and defined 



TV 

Co 



D{s] x,p) := —xf{s] x,p) 



j=l 3 



The radiated power is then 

-3 rt 



37r2 

Corollary 4. Let 



cj(t) := e-"-^d;cs/^cuoe"-^d;cs/^, 

where uq G J^i{QMX{H^ies)'^) > Banach space of trace class operators on 
QMX{H^Tes)'^ ! j'ar^za/ trace over the field states 

ujpit) := Tr^u;(t), 

then 

oopit) = e-<i/^^p(0)e"<p/^ + 
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where 

N 

2m,- 



^ 1 



N 



anc? denotes the space of trace class operators on M'^. 

Proof. The proof follows from the following three facts: 



1. the term of order e in equation (l9dl) is off-diagonal with respect to the 
Qm', 

2. the diagonal Hamiltonian defined in flU^ . is equal to p ® 1 + 
1 ® iff, so we have that 

tr^(e-'<V.^^e'<'/^) = e-'<i/^tr,^(cuo)e"<p/^ ; 

3. the following well known inequality, which holds for any Hilbert space 



□ 



Theorem 6. Let S be an observable for the particles, S G C{J^p), and 
uj G J^i{Pl^x{H'^)'^) Oj density matrix for a mixed dressed state with M free 
photons whose time evolution is defined by 

uj{t) := e-'*^'/^o^e'*^'/^ (102) 

then 

Tr^.((5® l^)^(t)) = Tr^^(5e-"<pTr^(a;)e"<p) + 

+ 0{e'/'mi - 6mo) + 0{eWiog{a{e)-^m + \t\')). 

(103) 

Proof. First of all we observe that, using proposition [3] and lemma [Tj we have 
Tr^ (^(5® l^)a;(t)^ = Tr^- (^(5 ® l^)e-'*^"*'/^cu,(,)■ 
■ e'*^"''/^^ + 0(a(£)i/2) , 
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where a;o-{e) G J^i(P|^x(-f^^''^)'^)- By the definition of the dressed Hamiho- 
nian and the cychcity of the trace we have then at the leading order 



Tr^( {S®l^)ujit) 



The transformed observable, using the definition of ^ and lemma [5l is given 
by 

® 1^) ^* =S®1^ + exCho)Ul, {S ® 1,^) + 
+ 0{e'\og{a-')). 

All the terms of order e in the previous expression are off-diagonal with 
respect to the Qms, and the same holds for the term of order e in ( |93l) . 
Therefore, they all vanish when we calculate the trace. Using point 2 and 3 
of last corollary we get then (11031) with 

instead of u. 

Using again lemma [T] and the fact that the terms of order e in the expan- 
sion of ^ are off-diagonal we can in the end replace by the identity and 
u^^^) by u. □ 

A The limit c ^ oo 

In this appendix we sketch the proof of theorem [H The reader can find 
additional discussions in |Dai| and ( |Sp| , chapter 17 and section 20.2). 



As remarked in the introduction, see equations (!9l)- (fT2l) . the limit c — > 
oo has the form of a weak coupling limit, in which the weak interaction 
is observed over the long time scale r = c^t. The corresponding physical 
interpretation is that the small system made up of the particles interacts 
with an environment (the quantized field) which is traced out to analyze the 
dynamics of the small system only. 
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The mathematical framework, as explained in |Dai| , whose notation is 
employed in this appendix, considers the Banach space 



^ := ^i(^), 

the space of trace class operators on with the corresponding trace norm. 

^ contains the convex subset of positive operators of trace one, which 
are the mixed states of the composite system particles plus field (density 
matrices). 

The Hamiltonian Hx defines on ^ what in semigroup theory is called an 
"implemented semigroup "0 ( [Xl] and references therein) via the usual formula 
for the Schrodinger evolution of the states 

1^/(^5) :=e-'*^-f?e+'*^^ . 

Contrary to what in general can happen for an implemented semigroup, Vj^ 
is strongly continuous on ^ ( ^Moj . theorem 2) and can therefore be written 
as 

where is the generator of V^'^, called the (total) Liouvillean. 
In the same way we define 

Ut{Q) :=e-'*^^^?e+"^^ t/^ =e-'*'-*, 

yA^^^^| ._ g-it(//f+A2/3Hp)^g+it(Hf+A2/3/fp)^ yA ^ g-it(Lf+A2/3Lp)_ 

and the Liouvilleans Li and L4/3 associated to the interaction Hamiltonians 
hi and /i4/3. 

If the Hamiltonian H which implements the group is a bounded operator, 
the corresponding Liouvillean L is also bounded and given by 

L{q) = [H,q], Mq^m . 

When H is unbounded a little more care is needed, because L will be also 
unbounded and defined only on a suitable dense domain (see, e. g., |PrTi] ) . 

One is interested in studying the dynamics of the particles in the limit 
A — > 0, given that at time t = the field is in the reference state cur G J^i(t^), 
which is invariant under the free dynamics Vt- Under these conditions a 
natural choice is = Qo? the projector on the Fock vacuum. Contrary to 
what happens for the limit e — *^ 0, it is not possible here to look at the case 
M 7^ 0, because there does not exist any density matrix in Qu^ which 
commutes with iff. 

"'^In this case an implemented group of isometries. 
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To implement these ideas mathematically one defines a projection on the 
particles states, 

where x is the characteristic function of the interval indicated and L is a 
fixed integer greater than 2. 

Lemma. P^'^ is a projection of norm 1. 

In the following we denote for simplicity Pq'^ simply by Pq. Moreover we 
put Pi := 1 — Pq and, given an operator A, we denote by A^^^^ := PiAPj. 
The dynamics of the particles are given then by 

W^:=P^V,^Po . 



t 1 

5 



We apply now a standard procedure to get an integral equation for 
called generalized master equation (see |Dai| or, in a more general context 
|AlLe] ■ section III.l). 

We denote by the diagonal evolution 

f/A _ g-it(U+AL(")+AV3L(ii))_ 

Applying Duhamel Formula and noting that L^''^'' = L^^g^ = 0, one gets 

V,' = Ut + \ f ds UU^r + + A^/^' f ds UUl'^;^ + Li)?)K' • 

Jo Jo 

Then 

< = Pof/' + X fds Uts^f'^PiK'Po + X'/' f ds f/^.Lf/g^PiK'^o, 
Jo Jo 

PiV,'Po = X f dsUlXi'^PoV^Po + X"^ f dsUlX^lJ^PoV^Po. 
Jo Jo 

Putting := PqU^, and introducing the new variables r = XH and 



and 



H ■— ^ O^t 

a = X'^u we get then 



Jo 

+ r t/aX^.(,_,)[iri,4/3(A, r-a) + K,/s,i{X, r - cr)]W^.,^+ 
Jo 

+ X^'' f ciaX^.(,_,)ir4/3,4/3(A, r - a)Wl 
Jo 



'A-2(7' 
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where 

K,,{\r):= dxX\Lf''^U^Lf'\ ^,j = 1,4/3 . 

Jo 

One can show that, when A 0, Kij{\, r) converges to 

/■oo 

so that 
where 

Spelhng out the terms in Ki^i and formulating the resuh in Jif instead 
of expressing it in ^ one gets theorem [TJ 
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